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Abstract 

A diffusion process for charge distributions in a phase space is ex- 
amined. The corresponding charge moves in a force field and under 
an action of a random field. There are the diffusion motions for coor- 
dinates and for momenta. In our model, an inner state of the charge 
is defined by a complex vector. The vector rotates with a great con- 
stant angular velocity with respect to the proper time of the charge. 
A state of the diffusion process is a (complex-valued) wave function 
on the phase space. As in quantum mechanics, we assume that, for 
the wave functions, the superposition principle holds. The diffusion 
process averages out vectors of inner states from different points of 
the phase space. A differential equation for this diffusion process is 
founded and examined. We demonstrate that the motion (described 
by this process) decomposes into a fast motion and a slow motion. The 
fast motion reduces an arbitrary wave function to a function from a 
subspace whose elements are parameterized by complex-valued func- 
tions of coordinates. The slow motion occurs in this subspace and 
it is described by the Schrodinger equation. The parameters of the 
suggested model are estimated. The duration of the fast motion is of 
order 10 -11 s. 



Introduction 



In [H E] , Maslov studied the problem of describing the motion of a charg 
distribution on a phase space under the action of the a random field on th 



charges. He assumed that the distribution depends on a hidden parameter, 
which varies under motion. In [I], it was proved under certain assump- 
tions (Theorem 3.2) that, if the initial charge distribution, which depends 
on coordinates and momenta, belongs to a certain subspace parameterized 
by complex-valued functions of coordinates, then the distribution does not 
leave this subspace, and the dynamics of such a system is described by the 
corresponding Schrodinger equation. The nature of the hidden parameter is 
inessential, but the dependence of the initial distribution on this parameter 
must be periodic with period of order h, where % is a small parameter (the 
Planck constant). 

Following pQ, we consider the diffusion of charge distribution, which de- 
pends on coordinates, momenta, and an angular parameter 9/h. We assume 
that the dynamics of the system with respect to coordinates and momenta 
is determined by a Hamiltonian and a heat diffusion process. 

The initial state can be any complex-valued function on the phase space 
(it is not required to belong to a certain function subspace, as in Theorem 3.2 
from |TJ). We show that, in the model under consideration, fast and slow 
motions can be distinguished. Under a fast motion, the systems starts in 
an arbitrary state and transfers into a state belonging to a subspace whose 
elements are parameterized by complex-valued functions of coordinates. A 
slow motion, as in Maslov's theorem, occurs in this subspace, and such a 
motion is described by the Schrodinger equation. 

In this paper, we also estimate the parameters of the suggested model. 

1 A description of the model 

We consider a mathematical model of a process whose state at each moment 
of time is determined by a wave function, that is, a complex-valued func- 
tion (p(x,p), where (x,p) G R 2n and n is the dimension of the configuration 
space. We assume that the change of states with time t in this procedure 
is determined by a certain (described below) Feynman path integral in the 
phase space. 

Consider the diffusion process on the phase space (x,p) G R 2n determined 
by a Fokker-Planck equation of the form 
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where f(x, p, t) is the probability density on the phase space at time t, H(x, p) 
is the Hamiltonian function, and a 2 and b 2 are the coefficients of the diffusions 
with respect to the coordinates and momenta, respectively. 

Diffusion process ([1]) determines in a standard way a probability measure 
V[x s , p s ]T>xT>p on the set of trajectories {x s ,p s } going from the point (x°,p°) 
at the initial moment s = and arriving at points (x,p) at time s = t. 

It is assumed that, if the wave function equals ip°(x,p) at the initial 
moment t = 0, then its value ip(x, p, t) at time t is determined by the following 
integral, which is an analogue of the Feynman path integral [3]: 

O, p) 

(f{x,p,t) =e abnt/h Jip°(x ,p°) J F[x s ,p s ]V[x s ,p s ]VxVp dx°dp°, (2) 

. t 

where F[x s ,p s ] = expj-^ J{H{x s , p s ) - x s p s )ds\, (3) 

o 

and the integral is with respect to the probability measure V[x s , p s ]T>xT>p 
on the set of trajectories {x s ,p s } defined above. The normalizing multiplier 
exp(abnt/h) of the integral is introduced for technical reasons. 

We emphasize that integral ([2]) under consideration differs from the Feyn- 
man integral [3] considered in quantum mechanics. In the case under consid- 
eration, the wave function depends on the coordinates and momenta, while 
in quantum mechanics, it depends only on the coordinates or only on the 
momenta. The quantum-mechanical Feynman integral is with respect to a 
"measure" defined by Feynman on the set of all trajectories in the configu- 
ration space. In the case under consideration, the integral is with respect to 
a probability measure on trajectories in the phase space. This probability 
measure is determined by the diffusion process (JTJ). 

As in quantum mechanics, we assume that, for wave functions, the su- 
perposition principle holds, and the probability density p(x,p) on the phase 
space that corresponds to a wave function <p(x,p) is defined by the standard 
formula 

p(x,p) = <p*(x,p)ip(x,p) = \ip(x,p)\ 2 . (4) 

Note that, like the Fokker-Planck equation, which is a continuous ap- 
proximate model (for sufficiently large times) of a more exact discrete model 
of a diffusion process, the integral representation (J2J) can be regarded as a 
continuous approximation of a more exact discrete model, which is described 
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by a complex Markov process. Such processes were studied by Maslov in|3] ■ 
In the case under consideration, the complex coefficients of the transition 
u(x l+1 , p i+1 ; x % ,p l \ At 1 ) from the point (x l ,p l ) to the point (x i+1 , p t+1 ) in time 
At 1 are given by 

u(x i+1 ,p i+1 ; x l , p i ; At 1 ) = e Zi K{x i+1 ,p i+1 ;x i ,p i ;At i ), 

nhwAt 1 i / \ 

where Z { = — ( H(x\ p l )At - (x i+l - x V ) , 

n n\ ) 

and K (x t+1 , p l+1 ; x 1 , p" 1 ; At 1 ) is the probability of the transition from (x\p l ) 
to , p l+l ) in time At 1 for a discrete diffusion process, which is described 
by the Fokker-Planck equation (TjTJ in the limit as At 1 tends to zero. 



2 Main resalts 

For integral (jSj), a differential equation is constructed by a standard scheme [3]. 

Theorem 1. The wave function determined by integral (TJ|) satisfies the 
differential equation 

dtp "fdHdp OH dp \ if TT "dH \ 

n ( d ipk\ 2 n d 2 abn 

where A ab <p = a 2 ^ z- ) <^ + & 2 ^ ^"T^ + (6) 

Removing the last term from Eq. (131) , we obtain a first-order partial differ- 
ential equation. This part of Eq. §5§ describes the deterministic component 
of the motion of complex vectors <p(x,p, t). According to the equation, under 
this motion, the application point of each vector moves along the classical 
trajectory determined by the Hamiltonian H(x,p), and the vector itself ro- 
tates with angular velocity to' = | \H — Ylk=i §^P^j a ^ eacn P om t of the 
trajectory. 

Note that, if the configuration space is three-dimensional and H = 
c\frri 2 c 2 +p 2 , then u'dt = v^m^JL — ^dr, where r is the proper time in 
the coordinate system related to a particle moving with momentum p. Thus, 
in this case, the vector, whose application point moves along a classical tra- 
jectory, rotates with constant angular velocity uj = mc 2 /h in the coordinate 
system related to this point. 
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Conversely, leaving only the last term of the form on the right-hand 
side of Eq. (J5J), we obtain the equation 



This equation describes the diffusion component of the motion of the vec- 
tors <p(x,p,t) in the phase space. Under this motion, the application points 
of vectors move in accordance with the classical homogeneous diffusion pro- 
cess. Under small displacements from a point (x, p) to a point (x + dx, p + 
dp), each vector is parallel-translated, and its length at time t is multi- 
plied by exp(abnt/h). The parallel translation of the vectors on the phase 
space is determined by the connection defined by L^x,dp)f{x,p) — <p(x,p) ~ 
— (i/h){p(x,p)pdq, where L( dXjdp )(p(x,p) is the translation of the vector <p(x,p) 
from the point (x,p) along the infinitesimal vector (dx,dp). In the special 
case where the configuration space is three-dimensional, such a connection 
on the phase space is caused by the synchronization of the moving clocks at 
the points of the phase space. The right-hand side of © is a self-adjoint 
operator. The eigenvalue problem for this operator reduces to the station- 
ary Schrodinger equation for harmonic oscillations by applying the Fourier 
transform with respect to coordinates. This implies that the eigenvalues of 
the operator of Eq. (J7J) are nonpositive, and the following theorem is valid. 

Theorem 2. Let (p(x,p, 0) be an arbitrary function whose Fourier trans- 
form with respect to p tends to zero as x — > oo. Then, the solution <p(x,p,t) 
to the diffusion equation |7|) tends to a stationary solution of the form 





<p(x,p) = lim <p(x,p,t) 

t^-OO 



1 



J^(y)x(x,y)e-^-^ h dy, (8) 



(27rh) n / 2 



where 





and 




(10) 



exponentially in time (with exponent 



abt/h). 
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Note that x 2 ( x >y) ls the probability density of the normal distribution 
in x with mathematical expectation y and variance ah /(2b). If the value of 
ah /(2b) is small, then the function x 2 ( x i y) is close to the delta-function of 
x-y. 

The composition of (Q and (jSJ) is a projector from the space of wave 
functions on the phase space onto some subspace. The elements of this 
subspace are parameterized by functions of the form ip(y), where y G R n , 
i e., by wave functions on the configuration space. 

Formulas © and (j3J) imply the following assertion. 

Theorem 3. Ifip(x) is a wave function on the configuration space, then 
the corresponding probability density on the phase space is determined by 

p(x, P ) = j^- J^(y)r(y')x(x,y)x(x,y')e-^-y'^ h dydy', (11) 

where x( x ,y) is defined by l[T(J\) . 

In contrast to the quasi-distributions defined by Wigner in [5], the density 
p(x,p) defined on the phase space by ffTTj) is always nonnegative. The algebra 
of observables determined by real functions on the phase space but averaged 
over probability distribution densities of the form (ITTj) was studied in [6]. 

The probability distribution density p(x) on the configuration space is 
expressed in terms of the density p(x,p) by the formula p(x) = f R3 p(x,p)dp. 
Using this expression and integrating (TTTT) with respect to p, we obtain the 
following assertion. 

Corollary 1. If ip(x) is a wave function on the configuration space, then 
the corresponding density p(x) on the configuration space is given by 

p(x) = f \4j(y)\ 2 X 2 (x,y)dy, (12) 

where x( x ,y) is defined by UP) . Thus, p(x) is obtained from \ip(x)\ 2 by 
smoothing (convolving) with respect to the density of the normal distribution 
with variance ah /(2b), and the accuracy of the determination of the coordi- 
nate is bounded by ~ ^Jah/(2b). 

As is known j7j, in quantum electrodynamics, the minimal error of mea- 
suring the coordinates of an electron in a system at rest is estimated as 
h/(mc), where m is the mass of the electron and c is the speed of light. 

In statistical physics (see, e.g., [8J, Chapter 7, Sections 4 and 9), under 
the assumption that diffusion is caused by thermal actions on electrons, the 
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diffusion coefficients in the coordinates and momenta are expressed in terms 
of the temperature T as a 2 = kT/{m r y) and b 2 = 'ykTm, where k is the 
Boltzmann constant, m is the mass of the electron, and 7 is the friction 
coefficient of the medium per unit mass. This implies a/b = (7m) -1 and 
ab = kT. Thus, in this case, the quantity a/b, which is contained in ffTUl) 
and determines the variance of smoothing in Corollary 1, does not depend on 
the temperature. Moreover, the duration t of the transition process, which 
is defined in Theorem 2, has the form t ~ h/(ab) = h/(kT) = T _1 • 7.638 • 
KT 12 s. 

Taking into account this estimate, we regard the quantity h/(ab) in Eq. (J5]) 
as a small parameter and assume that the coordinates and momenta little 
change during the transition time under the classical motion determined by 
the Hamiltonian H(x,p). 

Theorem 4. The motion described by Eq. |3]) asymptotically decomposes 
into a fast and a slow motion as h/(ab) — > 0. The fast motion reduces an ar- 
bitrary wave function (p(x,p, 0) to the form in time of order h/(ab). The 
wave functions of the form constitute a linear subspace. The elements 
of this subspace are parameterized by wave functions ip(y) depending only on 
the coordinates y G R n . The slow motion, which starts from a nonzero wave 
function belonging to this subspace, occurs in this subspace and is parame- 
terized by a wave function ip(y,t) depending on time. The function ip(y,t) 
satisfies the Schrodinger equation ihdip/dt = Hip, where 

£j , 1 /Yrr/ \ ^fOH ib8H\. ,A 

xxO, y)x(x, y')e^ (y ~ y ' )p ^{y\ t)dy'dxdp, 

and x( x ^y) is defined by /[TO]) . 

Theorem 5. If ^ is a small quantity and H(x,p) = + V(x), then, 
accurate to terms of order ah/b, the operator H has the form 

- h 2 /» a ! \ ah " d 2 V , 3nbh 

The first two terms in (JTBl give the standard Hamiltonian operator; the 
last term is a constant, and it can be neglected. The next-to-last term can 
be regarded as a perturbation of the Hamiltonian (because ah/b is small). 
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Assuming that the deviations in the spectrum of the hydrogen atom (the 
Lamb shift), which are observed in the Lamb-Retherford experiments [9], 
are caused by the next-to-last term in ([TBI , we can estimate the quantity 
a/b. Calculations by the standard method of perturbations, similar to those 
performed in [10], give the estimate a/b = 3.41 ■ 10 4 s/g. Thus, the stan- 
dard deviation for the normal distribution x 2 , over which smoothing in (TT21) 

is performed, has the form ^ ah/{2b) = 4.24 ■ 10~ 12 cm. This quantity is 
substantially smaller than the radius of the hydrogen atom and close to the 
Compton wave length h/(mc) = 3.86 • 10~ n cm of the electron. 

Thus, calculations show that the suggested model, which has the form of 
integral (|2J) or of differential equation ((Sj), sufficiently adequately describes 
physical processes in standard cases for the standard Hamiltonian. This 
model can also be applied to calculate processes with nonstandard Hamil- 
tonians or with Hamiltonians rapidly varying with time, which take place 
under sudden perturbations |Tl] or when the potential varies periodically 
with frequency of order ab/h; it would be interesting to compare the results 
with experimental data. 
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